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1. Introduction 

In the framework of the renormalization-group (RG) approach to critical phenomena, a quanti- 
tative description of many continuous phase transitions can be obtained by considering an effective 
Landau-Ginzburg-Wilson (LOW) field theory, containing up to fourth-order powers of the field 
components. The simplest example is the 0(/^)-symmetric theory, defined by the Lagrangian 
density 

^oiN) = \Y,{d^^if + yl^^l + i«(I<J>?)' (1-1) 

where is an A/^-component real field. These theories describe phase transitions characterized 
by the symmetry breaking 0{N)^0{N — 1). We mention the Ising universahty class for N = \ 
(which is relevant for the liquid-vapor transition in simple fluids, for the Curie transition in uniax- 
ial magnetic systems, etc.), the XY universality class for N = 2 (which describes the superfluid 
transition in "^He, transitions in magnets with easy-plane anisotropy, and in superconductors), the 
Heisenberg universality class for N = 3 (it describes the Curie transition in isotropic magnets). 
Moreover, the limit N ^0 describes the behavior of dilute homopolymers in a good solvent in the 
limit of large polymerization. See, e.g., Refs. |[T], ^ for recent reviews. 

Beside the transitions described by 0(A'^) models, there are also other physically interesting 
transitions described by more general Landau-Ginzburg-Wilson (LGW) O'* field theories, char- 
acterized by more complex symmetries and symmetry breaking patterns. The general LGW 
theory for an A/^-component field can be written as 

^ i ^ (■ ^- ijkl 

where the number of independent parameters r, and uij^i depends on the symmetry group of the 
theory. Here, we are only assuming a parity symmetry which forbids third-order terms. In the field- 
theoretical (FT) approach the RG flow is determined by a set of RG equations for the correlation 
functions of the order parameter. This approach has been applied to investigate the critical behavior 
of disorder and/or frustrated systems, magnets with anisotropy, spin and density wave models, 
competing orderings giving rise to multicritical behaviors, and also the finite-temperature transition 
in hadronic matter. 

The main issue discussed in this talk is the RG flow of general multi-parameter <I>^ field theo- 
ries, and its applications to the study of critical phenomena in statistical systems. In particular, we 
consider FT perturbative approaches based on expansions in powers of renormalized quartic cou- 
plings. An overview of FT results for physically interesting cases is presented and compared with 
experiments and other approaches such as lattice techniques. Finally, this RG approach is applied 
to the study of the nature of the finite-temperature transition in QCD with Nf light flavours. 

2. Renormalization-group theory of critical phenomena 

Critical phenomena are observed in many physical systems when they undergo second-order 
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Figure 1: Typical phase diagrams of ferromagnetic materials (left) and liquids (right) 



(continuous) transitions characterized by a nonanalytic behavior due to a diverging length scale. 
Classical examples are transitions in ferromagnetic materials and liquids, whose phase diagrams 
are sketched in Fig. |[ The first general framework to understand critical phenomena was proposed 
by Landau [^; it was essentially a mean-field approximation. A satisfactory understanding was 
later achieved by the Wilson RG theory [Q]. 

The main ideas to describe the critical behavior at a continuous transition are (i) the existence 
of an order parameter which effectively describes the critical modes; (ii) the scaling hypothesis: 
singularities arise from the long-range correlations of the order parameter, which develop a diverg- 
ing length scale; (iii) universality: the critical behavior is essentially determined by a few global 
properties, such as the space dimensionality, the nature and the symmetry of the order parameter, 
the symmetry breaking, the range of the effective interactions. The RG theory of critical phenom- 
ena [m5[ |6|] provides a general framework where these features naturally arise. It considers a RG 
flow in a Hamiltonian space. The critical behavior is associated with a fixed point (FP) of the RG 
flow where only a few perturbations are relevant. The corresponding positive eigenvalues of the 
linearized theory around the FP are related to the critical exponents v, T], etc... 

According to the RG theory, the singular part of the Gibbs free energy obeys a scaling law 

=^smg("l ,"2, • .■,Uk,...)= Z7"''=^sing(^^'"l>'-"2, • • • • • •) (2.1) 

where are the nonlinear scaling fields, which are analytic functions of the model parameters. In 
a standard continuous transition there are two relevant scaling fields with y, > 0, Ut and m/,, and an 
infinite set of irrelevant fields w, with yt < 0. The relevant fields may be identified with the reduced 
temperature t = T /Tc — I and the external (magnetic) field H, i.e. Ut ^ t and m/, ~ // for f,// ^ 0. 
Setting b^'' |Mf | = 1, we can write 

fsmg=^''' =^sing {Uh I Ur \ -^'1^ , Wi\u, j ->'/-^' ) (2.2) 

'Continuous transitions are generally characterized power-law behaviors. For example, defining the reduced tem- 
perature as ? = r/Tf — 1, in the disordered phase the magnetic susceptibility x and the correlation length £, diverge as 
X ~ t^^ and ^ t^^ respectively, the specific heat behaves as C// ~ kl the two-point function lV(^'(x) decays as 
1 /x''^^+'' at r = Tc, the magnetization vanishes as M ~ {—t)^ in the ordered phase, etc... In standard continuous transi- 
tions only two critical exponents are independent, because there are several scaling relations: y = (2 — T})v, a = 2 — dv, 
j8 = v(rf-2 + 77)/2. 
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Since w,|m, | ' ^ for f ^ 0, we can expand with respect to the ai^guments containing the 
irrelevant scahng fields, obtaining 



where = 1/v, j/, = (t/ + 2 — T] )/2, A,- = —yi/yt > 0. The first term of the r.h.s. of Eq. (23) is 
the universal asymptotic critical behavior, while the other terms give rise to nonuniversal scaling 
corrections. The above scaling relations can be easily extended to allow for finite-size systems. 
In the presence of other relevant perturbations beside t and H, one observes more complicated 
multicritical behaviors. 



3. Field-theoretical perturbative approach 

The RG theory provides the basis for the FT approaches to the study of critical phenomena. 
The critical behavior can be determined by the RG flow of a corresponding Euclidean quantum 
field theory (QFT). 

3.1 Critical behavior of statistical systems and RG flow of Euclidean QFTs 

Let us consider the Ising model defined on a (i-dimensional lattice: 

H = -J'£a,ay, a, = ±l, Z = ^ exp(-///r), (3.1) 

where the sum in the Hamiltonian is over the nearest-neighbor sites of the lattice. The critical 
behavior is due to the long-range modes, with / ^ 1. In order to describe these critical modes, 
one may perform blocking averages over sizes a <C /, and derive an effective Hamiltonian for the 
block average variables, which can be considered as real variables cp^ attached to the blocks. The 
simplest effective Hamiltonian of the block variables (px which preserves the Z2 symmetry may be 
written as 

H,s = a'-^Y.('P^+^^ - "P-O" + «a^£(<P.v' - v')' (3-2) 

This blocking procedure does not affect the long range modes at the scale /, and therefore it does 
not change the universality class of the transition. Then, we may consider the limit a ^ of H^ff. 

H^. = j d^x^{(p), j^((p) = l{d^(p)^ + ^r(p^ + ^u(p^ (3.3) 

where r — T — Tc. Again, this limit does not change the universality class. The corresponding 
partition function 

Z = J [d(p] exp[- J d''xJ^{(p)] (3.4) 

is the path integral of an Euclidean QFT with -Sf (<p) = J^{(p). Therefore, the critical behavior of 
the original Ising model is related to the behavior of the correlation functions of the QFT for 
a real one-component field in the massless limit. Analogous euristic arguments can be applied to 
other statistical systems and corresponding QFTs. 

Note that the way back provides a nonperturbative formulation of an Euclidean QFT, from the 
critical behavior of a statistical model. An important example is the Wilson lattice formulation of 
QCD defined in the critical (continuum) limit of a 4D statistical model. 
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3.2 Perturbative schemes in the field-theory approach 



A successful approach to the study of critical phenomena exploits the relation with the RG 
flow of a corresponding QFT, and therefore the typical techniques used in QFTs. We are interested 
in the critical behavior of the "bare" correlation functions r„(p;r,M,A), which correspond to the 
"physical" correlation functions of the statistical system, unlike high-energy physics where the 
physical correlation functions are the renormalized ones. 

In this section we consider the simplest case of the □(A'^) field theory ( |1.1D . For d <4, i.e. 
J = 3,2, the theory is super-renormalizable since the number of primitively divergent diagrams is 
finite. Such divergences are related to the necessity of performing an infinite renormalization of the 
parameter r appearing in the bare Lagrangian, see, e.g., the discussion in Ref. [^. This problem 
can be avoided by replacing the quadratic parameter r of the Lagrangian with the mass m (inverse 
second-moment correlation length) defined by 



m 



-2 



1 



5r(2)(;,2 



r(2)(o) dp^ 



(3.5) 



where rfj(p2) =r(2)(p2) 

5ah is the one-particle irreducible two-point function. Perturbation the- 
ory in terms of m and the bare quartic coupling u is finite in J < 4. In the massive zero-momentum 
(MZM) scheme, one considers a set of zero-momentum renormahzation conditions in the mas- 
sive (disordered) phase. 



T^^\p)=Z-'[m'+p^ + 0{p' 



r(4)(o) 



(3.6) 



where F^"^ are one-particle irreducible correlation functions, and g is the renormalized MZM cou- 
pling. Moreover, one defines = r(i'2)(0), where T^^-^\p) is the one-particle irreducible two- 
point function with an insertion of xhe critical limit m^O can be studied by writing Callan- 
Symanzik RG equations for the renormalized correlation functions f'"'' = Z^^^F^") , 

[2-T7^(^)]mM''"'(;,;0), 



T^:'\p) 



(3.7) 



see, for example, Ref. [0]. The RG functions 



^{g)=m—, 
dm 



dlnZ^, 
dlnm 



(3.8) 



can be computed as power series of g. For the 3D 0(A'^) <I>^ models they have been computed to six 
and seven loops respectively [[l0|]. When m — > the coupling g is driven toward an infrared-stable 
fixed point (FP), i.e. a zero g* of the j8-function, j8(g) = -(o{g* - g) + 0[{g* - gY] with x > 1. 
The Callan-Symanzik j3 -function fi{g) is nonanalytic at ^ = ^* [jlTl]. Using the RG equations, one 
can also identify Tj = ri^{g*) and 1/v = 2 - ri^{g*) + T],(^*). 



One may also consider an alternative perturbative scheme: the MS renormalization scheme 1 13], 
defined at T = T^, i.e. in the massless theory. This is based on the dimensional regularization and 
the subtraction of the l/e poles (e =4 — d) to obtain the renormalized correlation functions. One 
sets (£> = [Z0{g)Y^^(£>r, u = AdPL^Zg{g), where is a J-dependent constant, and g is the MS 
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3D Ising exponents 



a 



EXPT 



PFT 



liquid-vapour 
fluid mixtures 
uniaxial magnets 



6,7-/ MZM [|16p 
0{e^) exp [jT|] 



Lattice HT exp [|17|] 



MC [18] 



0.6297(4) 

0.6297(7) 

0.6300(17) 

0.6304(13) 

0.6290(25) 

0.63012(16) 

0.63020(12) 



0.111(1) 

0.111(2) 

0.110(5) 

0.109(4) 

0.113(7) 

0.1096(5) 

0.1094(4) 



0.042(6) 
0.038(3) 

0.034(3) 
0.036(5) 
0.0364(2) 
0.0368(2) 



0.324(2) 

0.327(3) 

0.325(2) 

0.326(1) 

0.326(3) 

0.3265(1) 

0.3267(1) 



Table 1: Estimates of the critical exponents of the 3D Ising universality class, from experiments (taken 
from the review resummation of the FT 6,7-loop calculations within the MZM scheme and of 0{e^) 
expansions, and from lattice techniques: 25th order high-temperature (HT) expansion and Monte Carlo (MC) 
simulations. 



renormalized quartic coupling. Moreover, one defines a mass renormalization constant Z, by re- 
quiring ZfP^''^) to be finite when expressed in terms of g. Then one derives the RG functions 
Pis) = l^dg/dfJ. and T70,r(g) = dlnZ^ t /dlnjj.. Analogously to the MZM scheme, the nontrivial 
zero g* of the jS -function is the stable FP, and the critical exponents are obtained by evaluating the 
RG functions T]^ , (g) atg =g*. Within this scheme, one can perform the so-called e=4 — d expan- 
sion [14], i.e. and expansion about d = 4. Alternatively, one can fix d after renormalization [p^]. 



for example d = 3, obtaining a fixed-dimension expansion in the MS coupling g. 

FT perturbative expansions are divergent. If we consider a quantity S{g) that has a per- 
turbative expansion S{g) » T^s^g'^, the large-order behavior of the coefficients is given by Sk ~ 
k\ {—a)^k!' [l + 0{k^^)\ , with a > for (i < 4. Thus, in order to obtain accurate results, an appro- 
priate resummation is required before evaluating the RG functions at the fixed-point value g* . This 
can be done by exploiting their Borel summability, which has been proved for the fixed-dimension 
expansion in J < 4, see e.g. Refs. [^ |]] and references therein, and has been conjectured for 
the e expansion. Moreover, one can also exploit knowledge of the high-order behavior of the ex- 
pansion, which is computed by semiclassical instanton calculations, see, e.g., Refs. [12, ^. Note 



that the results of this resummation are essentially nonpertubative, because it uses nonperturbative 
information. 

3.3 The 3D Ising and XY universality classes 

The effectiveness of the RG approach based on <I>^ QFTs can be appreciated by looking at the 
results for the 3D 0(A'^) models, and in particular for the Ising and XY universality classes. 

The Ising universality class corresponds to a theory with a real one-component field. It 
describes transitions in several physical systems, such as liquid-vapor systems, fluid mixtures, 
uniaxial magnets. Table [| shows selected results from experiments, field-theoretical and lattice 
computations. A more complete review of results can be found in Ref. [|l]]. FT results are quite 
precise and in good agreement with experiments and lattice computations, which provide the most 
precise theoretical estimates. 

The 3D XY universality class is characterized by a two-component order parameter and the 



symmetry U{\). It corresponds to the 0(A'^) model (1.1) with N = 2. An interesting represen- 
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3D XY exponents 



a 



EXPT 
PFT 

Lattice 



+He [19] 



6,1-1 MZM [jig 
0{£^) exp 0] 
MC+HT [|2 
MC[|l|] 



-0.0127(3) 

-0.011(4) 

-0.004(11) 

-0.0151(3)* 

-0.0151(9)* 



0.6709(1)* 

0.6703(15) 

0.6680(35) 

0.6717(1) 

0.6717(3) 



0.035(3) 
0.038(5) 
0.0381(2) 



Table 2: Results for the critical exponents of the 3D XY universality class, from experiments on "^He, 
resummation of the FT 6,7-loop calculations within the MZM scheme and of 0{e^) expansion, and from 
lattice techniques: a sinergy of HT expansions and MC simulations (MC+HT) and MC simulations (MC). 
Results marked by an asterisk are obtained by using the hyperscaling relation a = 2 — dv. 



tative of this universality class is the superfluid transition of ^He along the A -line T^{P), where 
the quantum amplitude of helium atoms is the order parameter. The superfluid transition of '^He 
provides an exceptional opportunity for a very accurate experimental test of the RG predictions. 
Moreover, experiments in a microgravity environment, for instance on the Space Shuttle [|^], can 
achieve a significant reduction of the gravity-induced broadening of the transition. Exploiting these 
favorable conditions, the specific heat of liquid helium was measured to within a few nK from the 
A -transition [|T9|]. In Table ^ we show selected results from experiments, field-theoretical and lat- 
tice computations. We note that there is a significant difference between the experimental result 
and the best theoretical estimates using lattice techniques. FT results are in good agreement with 
experiments and lattice computations, but they are not sufficiently precise to distinguish them. This 
discrepancy should not be necessarily considered as a failure of the RG theory, however it calls for 
further investigations. A proposal of a new space experiment has been presented in Ref. [|2^]. 



4. The RG flow in multi-parameter 4>'^ field theories and the 7] conjecture 

Beside transitions belonging to the 0(A'^) universality classes, many other physically inter- 
esting transitions are described by more general multi -parameter field theories, cf. Eq. (1.2), 



characterized by more complex symmetries [23, 0, |I|]. The number of independent parameters r, 
and Uijki depends on the symmetry group of the theory. An interesting class of models are those in 
which ^, is the unique quadratic polynomial invariant under the symmetry group of the theory, 
corresponding to the case all field components become critical simultaneously. This requires that 
all r; are equal, r, = r, and Uij^i must be such not to generate other quadratic invariant terms under 



RG transformations, for example, it must satisfy the trace condition [g5|] Y^iUuki °^ Ski- In these 
models, criticality is driven by tuning the single parameter r, which physically may correspond to 
the reduced temperature. More general LGW theories, which allow for the presence of indepen- 
dent quadratic parameters r, , must be considered to describe multicritical behaviors where there are 
independent correlation lengths that diverge simultaneously, which may arise from the competition 



of distinct types of ordering. Note that, like the simplest 0(A'^) models (1.1), all multi-parameter 
<I>^ field theories are expected to be trivial in four dimensions. 



7 



RGflow of multi-parameter field theories 



Ettore Vicari 



ordered 1 



0(n,) 



_flop line y disord 

phase 
ordered 2 \ 

0(n,) 



ordered 1 




ordered 2 



0(n,) 



disord 

phase 

0(112) 



T T 
Figure 2: Phase diagrams with muhicritical points where different transition lines meet. 

4.1 Physically interesting multi-parameter field theories 

Some physically interesting examples of multi-parameter <!>'* field theories are: 

• The O(M)0O{N) models, which is the most general O'* theory with aN xM matrix funda- 
mental field <t>ai {a= l,...,M,i= l,...,N), and the symmetry OiM)(^0(N), 

a,h i 

They are relevant for transitions in noncollinear frustrated magnets and in ^He. See Sec. ^ 

• The MN models with a real M xN matrix field Oq,-: 

(4.2) 



(4.1) 



i,ci ij.ab 



They describe transitions in randomly diluted M-component spins for N ^ 0, see below in 
Sec. ^ and magnets with cubic anisotropy for M = I, N = 2,3, see e.g. Refs. [ 23, ^ ]. 



The spin-density wave model has two complex A/^-component order-parameter fields: 

It is relevant for transitions where spin-density waves play an important role, as, for example. 



in high-Jc superconductors (cuprates), see e.g. Refs. [E7, 28]. 



Physically interesting examples of multicritical behavior arise from the competition of or- 
derings with symmetries 0(«i) and 0(«2)> at the multicritical point where the ti^ansition 
lines meet, as shown Fig. ^ The corresponding LGW <I>^ theory must have the symmetry 
0(ni)©0(n2) and two 0{ni) and 0(?i2) vector fields and ^2, 

The multicritical behavior is observed by tuning two relevant scaling fields, related to ri and 



^2 V-%- Multicritical behaviors of this type occur in several physical contexts: in anisotropic 
antiferromagnets, high-r^ superconductors, multicomponent polymer solutions, etc., see 



Refs. m, 30, 31,B3] 
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Figure 3: Example of RG flow in the coupling space of a two-parameter field theory. In particular, the 
figure sketches the RG flow in the renormalized coupling u-v plane of the 0{M)®0{N) <t>^ theory in the 
lai-ge-A^ limit, which has four FPs: the Gaussian (G), 0(M x A^) (H), chkal (C) and antichiral (A) FPs; the 
stable FP is the chiral one. 



4.2 Perturbative FT approach and RG flow 

The RG flow, and other interesting quantities as the critical exponents, can be determined by 
using perturbative approaches, extending the methods employed in the case of the □(A'^) models, 
see Sec. |3.2[ In the massive (disordered-phase) MZM scheme, one expands in powers of the MZM 
quartic coupUngs gijki , defined by 

Tf{p) = 8ijZ-^' [m^ + p^ + 0{p')\ , r|J(0) =mZ^-2^,,-,, (4.4) 

The massless (critical) MS scheme is based on a minimal subtraction procedure within the dimen- 
sional regularization, and can give rise to an s = 4 — d expansion, and also 3D expansions in the 
renormalized MS couplings gij^i by setting £ = 1 after renormalization. 

High-order computations have been performed for several LGW <I>^ theories, see e.g. Refs. |33, 



|26| , |34| , 1^, ^ 30, 37, 39, 38, 4^, 41, ^], to five or six loops, which requires the calculation of 
approximately 1000 Feynman diagrams. Again, the resummation of the series is essential. It can 
be done by exploiting Borel summability and calculation of the large-order behavior [p^]. This is 
achieved by extending the techniques employed for O(A^) models. The comparison of the results 
obtained from the analyses of the MZM and MS expansions provides nontrivial checks. 

The RG flow is determined by die FPs, which are common zeroes g?^; of the j8 -functions, 

PiM8abcd)=m^p^ (MZM), p.j,i(g^,^,)=^^pE (MS), (4.5) 
am ojX 

in the MZM and MS schemes respectively. A FP is stable if all eigenvalues of its stability matrix 
Sij = dpi/dgj\g=g* have positive real part. Fig. ^ shows an example of RG flow in the quartic- 
coupling plane of a two-quartic-parameter field theory. 

The existence of a stable FP implies that (i) physical systems with the given global proper- 
ties can undergo a continuous transition, (ii) the asymptotic behavior in continuous transitions is 
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controlled by the stable FP (apart from cases requiring further tunings). The absence of a stable 
FP predicts first-order transitions between the disordered and ordered phases of all systems. Note 
that, even in the presence of a stable FP, first-order transitions can be observed in systems that are 
outside the attraction domain of the stable FP. 

4.3 The T] conjecture 

Multi-parameter theories have usually several FPs. An interesting question is whether a 
physical quantity exists such that the comparison of its values at the FPs identifies the most stable 
FP. In 2D unitary QFT, such a quantity is the central charge c, which is related to the correlation 
function of the stress tensor 

(J{z,)T{z2)) = ^, (4.6) 

and to the behavior of the free-energy of a slab, which is /(L) \lx<=o = fooL — c'k/{6L) at T = T^. 
The c-theorem [ |42| ] implies that the stable FP in a 2D unitary QFT is the one with the least value 



of c. But, despite several attempts and some progress, see, for example, Refs. [|43[], no conclusive 
results on the extension of this theorem to higher dimensions has been obtained yet. 

Within general theories, the following conjecture has been put forward [Q]: In general 
unitary theories the infrared stable FP is the one that corresponds to the fastest decay of corre- 
lations. This corresponds to the FP with the largest value of the critical exponent 77 which charac- 
terizes the power-law decay of the two-point correlation function IV '^'(x) at criticality, 

w(^)Woc-^. (4.7) 

The exponent 77 is related to the RG dimension of the field, d,^ = {d — 2-\- r]) /2. Since in <i < 4 
general theories may have more than one stable FP with separate attraction domains, the T] 
conjecture should be then refined by comparing FP that are connected by RG trajectories starting 
from the Gaussian FP: among them the stable FP is the one with the largest value of t]. 

The 77 conjecture holds in the case of the 0(A'^)-symmetric O'* theory. For d <A, the Gaussian 
FP, for which T] = 0, is unstable against the non-trivial Wilson-Fisher FP for which T] > (the pos- 
itivity of T] in unitary theories follows rigorously from the spectral representation of the two-point 
function [^). It has been proven within the e = A — d expansion, i.e. close to d = A. It remains 
a conjecture at fixed dimension d <A, but its validity has been confirmed by several analytic and 
numerical results in lower dimensional cases. Several checks are reported in Ref. [p4|]. No coun- 
terexample has been found in all theories studied so far. This conjecture has also been extended 
to <I>^ theories describing multicritical behaviors, such as the 0(«i)©0(?i2) <I>^ theory (pT3|). In 
this situation, the exponent T] is replaced by a matrix and the conjecture applies to the trace of the 
matrix, i.e. TrT]. 

5. The 0{M)®0{N) 4>4 field theory 



In this section we discuss the RG flow of the 0{M)®0{N) theory ( |4~T| ). For M = 2 and 
vo > 0, the symmetry breaking is 0(2)(g)0(A'^)— 5-0(2)(g)0(A'^ — 2), which may be roughly simplified 
to 0{N)^0{N — 2) (although the 0(2) groups in the symmetry-breaking pattern are not the same). 
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These cases describe transitions in frustrated spin systems with noncollinear order [45], where 
frustration may arise either because of the special geometry of the lattice, or from the competition 
of different kinds of interactions. Typical examples of systems of the first type are stacked triangular 
antiferromagnets (STA's), for example CsMnBrs, CsVBrs, where magnetic ions are located at each 
site of a three-dimensional stacked triangular lattice. In STA's frustration gives rise to a further 
degeneracy of the ground state, which presents a chiral structure. This can be easily seen in the 
two-component XY antiferromagnetic model H = Y.{xy) ■ Sy defined on a triangular lattice. The 
minimal energy configurations have the chiral 120° structure shown in Fig. Another interesting 
case is given by M = 2,A'^ = 3, vq < 0, which is relevant for the superfluid transition of ^He, see 
e.g. Ref. @]. 

The nature of the transition in frustrated systems with noncollinear order has been longly 
debated. Experiments on several physical systems show continuous transitions, with exponents v = 
0.57(3), 0.54(3) for = 2 and v = 0.62(5) forA^ = 3(see e.g. Refs. [gl, [l|, |6[ |7|] for discussions 
of experimental results). The existence of these new chiral universality classes requires a stable 
FP in the 3-d RG flow of the corresponding 0{2)®0{N) <I>^ theory. This issue was investigated by 
using the e = A — d expansion. Calculations within the e = 4 — d expansion do not find any stable 
FP close to d = A, see e.g. Refs. |^, |], and references therein. The extension of this result 
to the physical dimensions d = 3 would predict first-order transitions for any system, in apparent 
contradiction with experiments. But the validity of such an extension to J = 3 is not guaranteed, 
because new FP's may appear going from d <A to d = ?>. Fixed-dimension d = 3 computations 
are required to check this possibility. 



The issue has been solved by high-order computations within 3D FT schemes [ |36[ , [49| , |47[ ] : 
6-loop and 5-loop in the MZM and MS schemes respectively. Their analyses provide a robust 
evidence of the presence of a stable FP, supporting the existence of new 3D chiral universality 
classes, which explain experiments. The left Fig. ^ shows the zeroes of the MS j3 -functions. The 
right Fig. || shows the RG trajectories in the u,v plane from the unstable Gaussian to the stable 
chiral FP (which can be obtained by solving the RG equation —Xdgi/dX = Pi{gj), with A G [0,°°) 
and gj{0) = 0, dgi/d?i\^^Q = Ui). Critical exponents V = 0.57(3) for = 2 and v = 0.55(3) for 
= 3 are in substantial agreement with experiments. The existence of these chiral universality 
class has been confirmed by MC simulations of lattice models [47, ^ ]. 

The RG flow of the 0{M)(^0{N) theory can be also studied analytically in the large-A^ limit 
keeping fixed M, for any 2 < d < A. In the \wge.-N limit one finds four FPs as shown in Fig. ^, 
the stable FP is the chiral one denoted by the letter C. The values of the critical exponent rj at 
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Figure 5: Some results for the RG flow of 3D 0{2)®0{N) models obtained by the analysis of the five- 
loop MS series The left figure shows the zeroes of the MS )3 functions. The right figure shows RG 
trajectories in the m,v plane from the unstable Gaussian and 0(2N) FPs to the stable chiral FP for various 
values of the ratio s = vq/mq of the bare quartic parameters. 



the FPs provide a further confirm of the Tj conjecture [44|, indeed setting rj = rjied /N + 0{l/N^) 
where ed is a constant depending on d, one finds [51] Tji = 0, r/i = 1/M, Tji = (M+ l)/2, and 
T]i = (M- l)(M + 2)/(2M) respectively for the Gaussian, 0{M xN), chiral ad antichiral FPs. AH 
numerical results at fixed M,N satisfy the tj conjecture. 

6. Ferromagnetic transitions in disordered spin systems 

The FT approach allows us to also successfully describe ferromagnetic transitions in disor- 
dered spin systems, which are of considerable theoretical and experimental interest. Such transi- 
tions are observed in spin systems with impurities, such as mixing of antiferromagnetic materials 
with non magnetic ones, for example Fe„Zni_„F2, Mn„Zni_„F2 (uniaxial), Fe^Er^, FevMnyZrj, 
(isotropic), and also "^He in porous materials. See e.g. Refs. [ p4l , | 



56] for experimental and 



theoretical reviews. These systems can be modeled by the lattice Hamiltonian 

Hp = -J Y,PxPv^x-Sy 



(6.1) 



where the sum is over nearest-neighbor sites, are M-component spin variables, and = 1,0 with 
probability p and I — p, respectively. The disorder is quenched: it mimicks the physical situation 
in which the relaxation time of the diffusion of impurities is much larger than other typical scales. 
This implies that the free energy F{p) <x InZ(p) must be averaged over the disorder. Accordingly, 
the expectation value of an observable must be computed by 



{P}) 



[dpnpmiPM) 



In the FT approach randomly-dilute spin models can be described by a field theory for an 
M-component field <I>; with and external random field coupled to the energy-density operator: 



(6.2) 
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Figure 6: The RG flow of the MN model for M = 1 and ^ in flie quartic-coupling plane. It shows the 
trajectories from the unstable Gaussian and Ising FPs to the stable RDIs FP for various values of the ratio 
s = Mo/^0 of the bare quartic parameters. Results obtained by the analysis of the six-loop MZM series [^3|]. 



where ^lf{x) is a spatially uncorrelated random field, with probability P(v^) ~ exp(— i//^^/4w). Then, 
using the replica trick, InZ = lim„^o(Z" — one can formally integrate over the disorder vari- 
ables, arriving at a translation invariant theory, which is the MN model ( t .2 ) with mq = — 6w < 0. 
The critical behavior of the original system is determined by the RG flow of the MN model in the 
nonunitary limit N ^ 0. Therefore, it can be determined by analyzing the high-order MZM and 
MS series for N = 0, which have been computed respectively to six loops [^, 35, 34] and to five 



loops |g3|,gO|. 

An interesting physical issue is whether the presence of impurities, and in general of quenched 



disorder coupled to the energy density, can change the critical behavior. General RG arguments [52, 
23] show that the asymptotic critical behavior remains unchanged if the specific-heat exponent a 
of the pure spin system is negative, which is the case of multicomponent 0(M)-symmetric spin 
systems, i.e. M > 1. On the other hand, a different critical behavior is expected in the case of 
Ising-like systems (M = 1), due to fact that otismg = 0.1096(5). This is confirmed by the results 
of the analyses of the high-order FT perturbative series, which show that the RG flow of the MN 
model in the limit N —s-O has a stable FP in the region m < when M = 1, see Fig. ^, implying the 
existence of a 3D randomly-dilute Ising (RDIs) universality class. 

Experiments, see, e.g., Refs. [|^, |l|] and references therein, confirm this scenario. The asymp- 
totic critical behavior remains unchanged for multicomponent systems. Table ^ reports results for 



Ising-like systems: from experiments |54], the analysis of the six-loop FT expansion in the MZM 



scheme [35], and recent Monte Carlo simulations [57]. The global agreement is very good 



It is worth mentioning that the RDIs universality class also describes ferromagnetic transitions 
in the presence of frustration, when frustration is not too large. This is for exmaple found in the 3D 



±J Ising model []58|], defined on a simple cubic lattice by the Hamiltonian 



H±J = - ^ JxyOxOy, ajc = ±\, 



(6.3) 



where Jxy, = ±7 with probability P{Jxy) = p5{Jxy —/) + (! — p)5{Jxy + J)- This is a simplified 
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RDIs exponents 



EXPT 0] 

PFT 6-1 MZM @ 

MC[g7|l 



0.69(1) 

0.678(10) 

0.683(2) 



0.359(9) 
0.349(5) 
0.354(1) 



Table 3: Results for the critical exponents of the 3D RDIs universality class. 



model [p9|| for disordered and frustrated spin systems showing glass behavior in some region of 



their phase diagram. Unlike model (5.1), the ±/ Ising model is frustrated for any p. Neverthless, 



the paramagnetic-ferromagnetic transition line extending for \ > p > pi^ = 0.76820(4) belongs 



to the RDIs universality class, i.e., frustration turns out to be irrelevant at this transition []58[]. For 
p < Pn the low-temperature phase is glassy with vanishing magnetization, thus the critical behavior 
at the transition belongs to a different Ising-glass universality class. 

7. The finite-temperature transition in hadronic matter 

The thermodynamics of Quantum Chromodynamics (QCD) is characterized by a transition 
at r ~ 200 Mev from a low-T hadronic phase, in which chiral symmetry is broken, to a high- 
T phase with deconfined quarks and gluons (quark-gluon plasma), in which chiral symmetry is 
restored [^]. Our understanding of the finite-r phase transition is essentially based on the relevant 
symmetry and symmetry-breaking pattern. In the presence of Nf light quarks the relevant symmetry 
is the chiral symmetry SU(A'y )l (g) SU(A^/)r. At T = this symmetry is spontaneously broken to 
SU(A'y)y with a nonzero quark condensate (yy)- The finite-r transition is related to the restoring 
of the chiral symmetry. It is therefore characterized by the simmetry breaking 

SU(A^/)L®SU(A^/)R^SU(A^/)y. (7.1) 

If the axial U(1)a symmetry is effectively restored at T^, the expected symmetry breaking becomes 

UiNf)L®\J{Nf)R^\J{Nf)v. (7.2) 

A suppression of the anomaly effects at Tc is however unlikely in QCD. Semiclassical calculations 
in the high-temperature phase [^] show that instantons are exponentially suppressed for T » 
Tc, implying a suppression of the anomaly effects in the high-temperature limit. Some lattice 
studies [ |62| ] suggest a significant reduction of the effective U(1)a symmetry breaking ai^ound Tc, 
but not a complete suppression. Since the anomaly, d^J^ °^ ^Q, gets suppressed in the large-A'^c 



limit, the symmetry-breaking pattern ( |7.2[ ) may be relevant in the large-A^c limit. 

Other interesting QCD-like theories are SU(A'^c) gauge theories with Nf Dirac fermions in the 
adjoint representation (aQCD). They are asymptotically free only for Nf < 11/4, thus only the cases 
Nf = 1,2 are interesting. Unlike QCD, aQCD is also invariant under global Z/v, transformations 
related to the center of the gauge group SU(A^c)> as in pure SU(A^c) gauge theories. There are two 
well-defined order parameters in the light-quark regime, related to the confining and chiral modes, 
i.e. the Polyakov loop and the quark condensate. Therefore, one generally expects two transitions: 
a deconfinement transition at T^i associated with the breaking of the Zyy^ symmetry, and a chiral 
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U(1)a anomaly 


suppressed anomaly at T^. 


QCD 


SU(iVy)L ® SUiNf)R ^ SUiNf)v 


U{Nf)L0U{Nf)R^UiNf)v 


Nf = \ 


crossover or first order 


0(2) or first order 


Nf = 2 


0(4) or first order 


U(2)l®U(2)k/U(2)v or first order 


Nf > 3 


first order 


first order 


aQCD 


SU(2A^j) ^ S0(2A^j) 


U{2Nf) 0{2Nf) 


Nf = l 


0(3) or first order 


U(2)/0(2) or first order 


Nf = 2 


SU(4)/SO(4) or first order 


first order 



Table 4: Summary of the RG predictions. We report the possible types of transition for each case, indicating 
the universaUty class when the transition can also be continuous. 



transition at Tc in which chiral symmetry is restored. In aQCD with Nf massless flavors the chiral- 
symmetry group extends to [^3|] S\J{2Nf). At T = this symmetry is expected to spontaneously 
break to S0{2Nf), due to quark condensation. Therefore the symmetry breaking at the finite-r 
chiral transition is 

SU(2A^/) ^ S0(2A^/) (7.3) 

with a symmetric 2Nf x 2Nf complex matrix as order parameter related to the bilinear quark 
condensate. If the axial U(1)a symmetry is restored at the symmetry-breaking pattern is 
U(2A^/) 0{2Nf). MC simulations for A^c = 3 and A/y = 2 ||4|, |5]] show that the deconfine- 
ment transition at is first order, while the chiral transition appears continuous. The ratio between 
the two critical temperatures turns out to be quite large: T^/Td ~ 8. 

In order to study the nature of the finite-r chiral transition in QCD and aQCD, one can exploit 
universality and renormalization-group (RG) arguments. [66, 37, SI, 67] 



(i) Let us first assume that the phase transition at is continuous for vanishing quark masses. 
In this case the length scale of the critical modes diverges approaching T^, becoming even- 
tually much larger than 1 /Tc, which is the size of the euclidean "temporal" dimension at T^. 
Therefore, the asymptotic critical behavior must be associated with a 3D universality class 
with the same symmetry breaking. The order parameter must be an Nf x Nf complex-matrix 
field ^ij, related to the bilinear quark operators \j/uYRj- 

(ii) The existence of such a 3D universality class can be investigated by considering the most 
general LGW theory compatible with the given symmetry breaking, which describes the 
critical modes at Tc. Neglecting the U(1)a anomaly, it is given by 

^u(;v) =Tr(5^<J>^)(5M<J>) + '-Tr<I>^<I' + J (Tr<I>^<J>)' + jTr(<I>1"<D)'. (7.4) 

If is a generic N xN complex matrix, the symmetry is \J{N)l®\J{N)r, which breaks to 
\J{N)v if vo > 0, thus providing the LGW theory relevant for QCD with Nf = N. If <I>,-y is also 
symmetric, the global symmetry is U(A'^), which breaks to 0(A'^) if vq > 0, which is the case 
relevant for aQCD with Nf = N /2. The reduction of the symmetry to S\J{Nf)L^SU{Nf)R 
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Stable FP 

zeroes of p u 
zeroes of (3 v 




Gaussian 0(n) 
^ unstable 

Figure 7: Zeroes of the /3 -functions associated with the quartic couplings of the Lagrangian ( [7.4| ) for = 2 
(fromRef 

for QCD [SU(2A'y ) for aQCD], due to the axial anomaly, is achieved by adding determinant 
terms, such as 

^su(iv) = ^u(/v) + wo (det<D^ + detO) . (7.5) 

Nonvanishing quark masses can be accounted for by adding an external-field term Tr (//O + h.c. 

(iii) The critical behavior at a continuous transition is determined by the FPs of the RG flow: the 
absence of a stable FP generally imphes first-order transitions. Therefore, a necessary con- 
dition of consistency with the initial hypothesis (i) of a continuous transition is the existence 
of stable FP in the corresponding LGW <I>^ theory. If no stable FP exists, the finite- r chiral 
transition of QCD (aQCD) is predicted to be first order. If a stable FP exists, the transition 
can be continuous, and its critical behavior is determined by the FP; but this does not exclude 
a first-order transition if the system is outside the attraction domain of the stable FP. 

The above RG arguments show that the nature of the finite-J transition in QCD and aQCD 
can be investigated by studying the RG flow of the corresponding 3D LGW field theories. 
RG studies based on high-order perturbative calculations in the MZM and MS are reported in 
Refs. [ pTj , PI ^ ]. Table ^ presents a summary of the predictions obtained by these RG analyses 
for the finite-r chiral transitions in QCD and aQCD with Nf massless quarks. 



The case relevant to Nf = \ QCD, i.e. the Lagrangian (7.4) with = 1, reduces to the 0(2) 
symmetric theory, corresponding to the 3D XY universality class. The determinant term related 
to the axial anomaly, cf. Eq. ( |7.5| ), plays the role of an external field, thus no continuous transition 
is expected, but a crossover. 

In the case relevant for Nf = 2 QCD with suppressed U(1)a anomaly, i.e. the O'* theory 
( fO| ) with N = 2, the analyses of both MZM and 3D MS schemes provide a robust evidence of a 
stable FP, see Fig. ^ A corresponding 3-D U(2)(g)U(2)/U(2) universality class exists, with critical 
exponents v ^ 0.7 and rj f» 0.1. No stable FP is found close to (i = 4 by one-loop £-expansion 
calculations [56], but, as already remarked in Sec. ||, the extension to (i = 3 of £-expansion results 



may fail. A stable FP is also found when the field <I> is symmetric, which is relevant for aQCD 
with one flavor, thus showing the existence of a universality class characterized by the simmetry 
breaking U(2)^0(2). 
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Figure 8: The scaling function Eij) for the 0(4) universality class, from Ref. [69|. The different lines 
represent different approximations. 



In two-flavor QCD, taking into account the U(l)^ anomaly, the symmetry breaking ( |7.1[ ) be- 
comes equivalent to the one of the 0(4) vector universality class, i.e. 0(4)— >0(3). The symmetry 
breaking ( |7.3| ) ofNf = \ aQCD is instead equivalent to the one of the 0(3) vector universality class. 
This already suggests that, if the transition in the chiral limit of Nf = 2 QCD is continuous, it must 
show the same asymptotic behavior of the 3D 0(4) universality class (0(3) in the case of Nf = I 
aQCD). This implies that the critical equation of state of the 3D 0(4) universality class, 

MocH\H\^^-^^^^E{y), yoct\H\-^/^^+^\ (7.6) 

provides a scaling relation between the quark condensate (iffxif)' the quark mass nif, which 
correspond respectively to the magnetization M and the external field H. The critical equation of 
state of the 3D 0(4) universality class has been accurately determined in the 3D 0(4) vector model: 
the critical exponents are | |6^ ] d = 4.789(6), j8 = 0.3882(10), and the universal scaling function 
E{y) is shown in Fig. ||. 

Actually, the LGW O'* theory corresponding to Nf = 2 QCD is quite complicated [|^, 67]: 



=^su(2) = Tr(5^<I>' )(5^<I>) + rTr<I>t<I>+ (TrO+O)' + ^Tr (Ol-O)' + 
+W0 (detO"'' + detO) + ^ (Tr<I>"f<I>) (detO"^ + detO) + ^ [{det<^^f + (detO)^ 



(7.7) 



where wo,xo,yo ^ g and g parametrizes the effective breaking of the U(1)a symmetry. If the 
anomaly is suppressed (g = 0), then wq = xq = yo = 0. ^su(2) contains two quadratic (mass) 
terms, therefore it describes several transition lines in the T-g plane, which meet at a multicritical 
point for g = 0. In the case of QCD the multicritical behavior is controlled by the U(2)/,(8)U(2)s 
symmetric theory. Possible phase diagrams in the T-g plane are shown in Fig. ^ When g ^0 the 
transition may be first order or continuous in the 0(4) universality class. Actually, we may also 
have a mean-field behavior (apart from logarithms) for particular values of g, see Fig. ^ If [.g | is 
small (a partial suppression of anomaly effects around is suggested by MC simulations [|6^]) 
we may observe crossover effects controlled by the U(2)(g)U(2) multicritical point at g = 0: if the 
transition is continuous at ^ = then the free-energy should behave as ^sing ~ fis^ '^)^ where 
t oc T — Tc{g = 0), V f« 0.7, 1.5. If the transition is first order at g = 0, then it is expected to 
remain first order for small g. A similar scenario applies also to Nf = I aQCD. 
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Figure 9: Possible phase diagrams in the T-g plane for the LGW theory ( [7.7[ ) describing the transition of 
Nf = 2 QCD, in the case the transition at the multicritical point, i.e. for g = 0, is continuous (left) or first 
order (right). Thick black lines indicate first-order transitions. At their end points, thus for particular values 
of g, the transition should be of mean-field type (apart from logarithms). 



No stable FPs are found for > 2 in the LGW theory ( |7.4D . Thus, neglecting the anomaly, 
transitions are always first order when Nf > 2 for QCD and Nf > I for aQCD. In most cases 
this result does not change if we take into account the axial anomaly, cf. Eq. (7^). The only 
exception is the case related to Nf = 2 aQCD, where a stable FP is found [41], corresponding to a 
3D SU(4)/SO(4) universality class, with critical exponents v « 1.1 and rj « 0.2. 

In nature quarks are not massless, although some of them, the quarks u and d, are light. The 
physically interesting case is QCD with Nf = 2 light quarks and four heavier quarks (in particular 
the quark s with m., ^ 100 Mev). Therefore, it is important to consider the effects of the quark 
masses in the above transition scenarios. According to the above RG arguments, if the transition is 
continuous in the chiral limit then an analytic crossover is expected for nonzero values of the quark 
masses mf, because the quark masses act as external fields in the corresponding LGW <I>^ theories. 
On the other hand, a first-order transition is generally robust against perturbations, and therefore it 
is expected to persist for mf > 0, up to an Ising end point. Actually, the presence of the massive 
quark s makes the above scenario more complicated, because the nature of the transition may be 
sensitive to its mass nis. Since the transition is expected to be first order in the chiral limit ofNf = 3 
degenerate quarks, we also expect that the first-order transition persists for sufficiently small value 
of On the other hand, if the transition is continuous in the limit oo corresponding to Nf = 2 
degenerate quarks, then there must be a tricritical point at = m* (where the critical behavior is 
mean field apart logarithms) separating the first-transition line from the 0(4) critical line. 

The nature of the transition in QCD can be investigated by lattice MC simulations. For nif >0 
around their physical values, MC simulations show that the low-J hadronic and high-J quark- 
gluon plasma regimes are not separated by a phase transition, but by an analytic crossover where 
the thermodynamic quantities change rapidly in a relatively narrow temperature interval, see, e.g., 
Refs. [^, 71, ^ 0]. Neverthless, the nature of the transition in the chiral limit is still of 
interest. Since the physical masses of the lightest quarks u and d are very small, some scaling 
relations may still be valid at the physical values of the quark masses, such as, for example, the 
0(4) relation (7^) between the quark condensate and masses in the case the phase transition in the 
chiral limit is continuous and belongs to the 0(4) universality class. 

The numerical investigation of the transition in the chiral limit is a hard task because it must 
be studied in the infinite-volume limit {V —>■ o°), in the continuum limit (N °° where N is the 
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number of lattice spacings along the Euclidean time direction), and massless limit (my 0). A 
robust control of the scaling corrections related to the continuum limit is essential. One cannot even 
exclude that in some lattice formulations the nature of the transition changes when approaching the 
continuum limit. Moreover, in the case of a continuous transition, the scaling corrections in the 



r.h.s. of Eq. ( |2.3[ ) may hide the universal asymptotic behavior if MC simulations are not done 
sufficiently close to the critical point (for example, in the case of the 0(4) universahty class, the 
leading scaling-correction exponent is not large, i.e. A ~ 0.5). 

In the case of Nf = 2 light quarks, the RG prediction for the chiral transition is that it is 
first order or a continuous transition in the 3D 0(4) universality class. Many studies based on 



MC simulations of different lattice formulations of QCD have been performed, see e.g. Refs. [75 



[76| , [77| , [78| , 79, 80, 81], but this issue is still controversial. Some MC results favor a continuous 
transition. However, the results have not been sufficient to settle its 0(4) nature yet. There are 
also results indicating a first-order transition. Unlike Nf = 2 QCD, the transition scenario appears 



settled for Nf >3: MC simulations [|77|, |82|, ^ |4|, |73|] show first-order transitions, in agreement 



with the RG predictions. Finally, in the case of Nf = 2 aQCD the available MC results [ |64| , |65| ] 
favor a continuous transitions. But they are not yet sufficiently precise to check the critical behavior 
of the 3D SU(4)/SO(4) universality class. 
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